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Abstract 

The distance dependence of the probability of observing two pho- 
tons in a waveguide is investigated. The Glauber correlation functions 
of the entangled photons in waveguides are considered and the spatial 
and temporal dependence of the correlation functions is evaluated. 
We derive upper bounds to the distance dependence of the probabil- 
ity of observing two photons. These inequalities should be possible to 
observe in experiments. 

1 Introduction 

The transmission of light in waveguides, in particular its quantum properties, 
is a topic of great interest in optics. Investigations of quantum correlations 
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and entanglement among photons have been in focus of in the foundations 
of quantum theory and its applications to quantum information science and 
metrology. 

With the emergence of quantum communication links over long distances 
[H 121 El H] , there is a need for a detailed study of the dependence of correla- 
tions and entanglement among photons on distance. 

Since direct interactions between photons in free space are extremely 
weak, generation of correlated photons generally requires nonlinear media 
such as the parametric down conversion. Recently, studies of two-photon 
scattering from a two-level system inside a one-dimensional waveguide have 
reported various features of photon correlation [21 El ?, E]- In particular a 
formal scattering theory to study multi-photon transport in a waveguide was 
employed [Sj. We also mention a general model describing nonlinear effects 
in propagation to and from the system in the quantum state, based on one- 
dimensinal model of field-atom interaction [9], see even [ID] . It describes 
spatiotemporal quantum coherence for the case of spontaneous emission from 
a single excited atom. In [5] this model was applied to the two-photon input 
wave packets. This field of research, spatiotemporal behavior of correlations 
of two photons propagating in nonlinear media, is closely related to studies 
on nonlinear response of a single atom to an input of two photons, e.g., from 
a single photon source [IT] . This response can be observed in the correlations 
between the two output photons. Here it is also very important to understand 
spatiotemporal dependence of correlations. 

There are thus sound reasons to study how correlation functions of entan- 
gled photons in hollow waveguides behave in space and time. The physical 
mechanism to model is dispersion that spreads pulses in space and time caus- 
ing attenuation with distance. The spreading limits also the bit rate for a 
given waveguide length because of mixing of pulses. Repeaters can be used 
at some length intervals, which cause higher costs and problems with pre- 
serving the quantum state through the repeater. Information on space and 
time properties of the correlation functions is thus of engineering interest. 
Such information is provided in the current paper by modelling the effect of 
modal dispersion. 

The study is also a preparation for more elaborate models, including ma- 
terial dispersion in hollow waveguides and fibres, cf. e.g. [12], [13]. To 
describe the situation, a brief background is presented [H] on optic fibres. 
For a fibre, it is possible to purify the material to the extent that losses from 
scattering from the impurities can be neglected in some wavelength bands. 
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Examples of such bands are the 1.3 jim. and 1.55 jum bands. The dispersion 
can be reduced to a negligible level for a certain wavelength within such a 
band, and a low dispersion can then be reached for a narrow wavelength 
band. Such a reduction of dispersion can be reached [Hj with a dispersion 
compensator which is a special fibre with tuned length, after the transmis- 
sion fibre, having the opposite dispersion to the transmission fibre. Another 
method [H] is dispersion shift, i.e., a choice of the transversal fibre dimen- 
sions. The chromatic dispersion, defined as the combination of modal and 
material dispersion, can be made to vanish for both methods at the design 
wavelength. 

It is worth mentioning the existence of repeaters for classical optic fibres 
using optical rather than previously used electrical methods [T5] and fibre 
switches that preserve the quantum state of the photon [16J . 

We shall study the asymptotic behaviour of the Glauber correlation func- 
tions for the entangled states of two photons in waveguides and show their 
vanishing for large distances. We estimate the rate of decrease of correlations 
and present upper bounds for the correlation functions. 

To estimate the correlation function we shall use some results on the 
properties of solutions of the (l+l)-dimensional Klein-Gordon equation anal- 
ogously used in the Haag-Ruelle scattering theory. 

We prove that the probability density P(zi, ti, z 2 , t 2 ) observing one pho- 
ton at point Z\ along the waveguide at time t\ and another photon at point 
z 2 at time t 2 satisfies the following inequality 

for some to and all zi,t\, z 2 ,t 2 . 

The bound (CQ) and the bounds ( 1261) . (|27j) (see below) should be possible 
to observe in experiments. 

Note that the decreasing of the correlations for entangled states in empty 
space was found in [?], see also discussion in [?], [?]. In this paper we have 
considered the waveguides and found the universal bound for the correlations, 
see also [?] for preliminary considerations. 

2 The photon probability density 

Let Ej(r,t) be the j-th component (j = 1,2,3) of the electric field operator 
at the space time point r, t. The operator can be written as the sum of the 
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positive and negative frequency parts: Ej(r,t) = Ej + \r,t) +Ej(r,t). The 
probabilities of photo detection are given by Glauber's formulas, [17] . In par- 
ticular, the probability that a state tp of the radiation field will lead to the 
detection at time t of a photon with the polarization along the direction j by 
a detector atom placed at point r is proportional to the first-order correlation 
function P^,{r,t,j) = (ip\E i j ~ ) (r,t)E!- + \r,t)\^). The joint probability of ob- 
serving one photo ionization with polarization j\ at point r\ at time between 
t\ and t\ + dti and another one with polarization j 2 at point r 2 between t 2 
and t 2 + dt 2 with t\ < t 2 is proportional to P^,(r 1 ,t 1 , j\, r 2 ,t 2 , j2)dt\dt 2 , where 
the second-order correlation function is 

P^r u t 1 Jr,r 2 ,t 2 J 2 ) = {^;\r 1 ,t 1 )E^\r 2 ,t 2 )E^ 

(2) 

The waveguide is a hollow conducting cylindrical tube r C R 3 along 
the z axis with boundary surface S and with a cross section Q with the 
bounding curve dQ in the xy-plane. It will be assumed that the walls have 
infinite conductivity. Appropriate boundary conditions are posed: E t \g = 
0, H n \s = 0, where E t is the component of electric field E tangential to the 
boundary of the waveguide and H n is the component of magnetic field H 
normal to the boundary. 

It is well known that in the interior of the waveguide the solutions of the 
Maxwell equations without sources can be divided into two sets of solutions, 
the so called TM modes with H z = and the TE modes with E z = [191120]. 
A general solution is a linear combination of the TM and TE modes. 

The general solution of the Maxwell equations in the waveguide can be 
written as follows [18J. Let <p n „(z,t) be any function satisfying the Klein- 
Gordon equation 

Q 2 Q 2 

' +m 2 nu )v nu (z,t) = 0. (3) 



dt 2 dz 2 

Here n = 1, 2, ... and v = TM or TE. We define for v = TM 



d 

E nu (r, t) = e nv (x, y)m nu —<f ni/ (z, t) (4) 

+n z m nu v n (x, y)ip nu (z, t), 

d 

H nv (r,t) = -h w (x,y)m nu —(p ni/ (z,t) 
where v n is the solution of the eigenvalue problem 
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d 2 d 2 

+ TT^ + m n) V n(%, U) = 0, (x, y) E tt, (5) 



dx 2 dy 2 

v n \an = 

with the properties 



v n (x,y)v n >(x,y)dxdy = 5 nn >, (6) 
^2v n (x,y)v n (x',y') = 5(x - x')5(y - y'), (7) 

n 

where w? n > are the eigenvalues. 
It is defined for v = TM 

e nv (x, y) = V T v n (x, y), (8) 

h nv (x,y) = n z x V T v n (x,y) = n z x e nu (x,y). 

E nu (r,t) and H nu (r,t) are defined in an analogous manner [?] for v = TE. 

The general solution of the Maxwell equations in the waveguide can now 
be written in the form 

E(r,t) = y2 E nu(r,t),H(r,t) = x [^Y^H nv {r,t). (9) 

V A*o — 

We write the solution of the Klein-Gordon Eq. ([3]) in the form 

/Hh 
n> m «(k)e iUnik)t - ikz (10) 

+a n {k)e~ iuln{ - k)t+ikz ), 

where oo n (k) = \Jk 2 + m^, and quantize it by taking a n (k),a^(k) as the 
annihilation and creation operators. 

Now the quantum electromagnetic field in the waveguide is reduced to a 
set of massive (l+l)-dimensional Klein-Gordon fields. Let us consider one 
of the modes. We define a one particle state 

IV*> = / 9(k)a[\0), (11) 



where |0) is the Fock vacuum. The probability density to detect the photon 
at the point z along the waveguide at time t is proportional to 

P(z,t) = {^-\z,t)cpM(z,t)\^), (12) 

where 

and Uk = y/k 2 + m 2 , m > 0. 

The expression (fT2l) can be written as 

P(z,t) = \A(z,t)\ 2 , (14) 

where 

/_— iui^t+ikz 
dkg(k) 2 ^__ . (15) 

The function A(z,t) is a solution of the Klein-Gordon equation. Let us 
consider the question of how the solution behaves in the limit of large t in a 
frame of reference moving with a constant velocity v < V < 1. In this frame, 
corresponding to the time x(l — v)/v after the wave front at t — x, the field 
is given by 

/ e -it(u k -kv) 
dkg(k) 2 ^__ . (16) 

For sufficiently smooth function g(k) one can use the stationary phase method 
to get 

p -it(u) ko -k v)-iTr/4: / Wjr 1 

Here fco is the solution of the equation to (k) = v, i.e. fco = mv/\ / l — v 2 and 
one has u (ko) > 0. For the probability P(z, t) we obtain 

P(vt,t)= 1 - A lgik °)jl, +0(-L). (18) 

More elaborated results on asymptotic expansions of the solution to the 
Klein-Gordon equation are given by Hormander 
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3 The entangled photon correlation function 
with distance dependence 

Now we define a two particle entangled state (biphoton) 

\1>) = j /(A^ 2 )«|0>, (19) 

where |0) is the Fock vacuum and f(h, k 2 ) is the two-photon wave function 
which is a symmetric function, f(h,k 2 ) = f(k 2 ,h) because we deal with 
bosons. 

The probability to detect one particle at the point z\ along the waveguide 
at time t\ and another particle at the space point z 2 at time t 2 is proportional 
to 

P{z 1 ,t 1 ,z 2 ,t 2 ) = (^-\z u t 1 )^-\z 2 ,t 2 ) V (+ \z 2 ,t 2 )^ + \z 1 ,t 1 )\^j). (20) 
The expression for P(z 1 ,t 1 , z 2 ,t 2 ) (1201) can be written as 

P(zi,ti,Z2,t 2 ) = \A{z x M,Z2,t 2 )\ 2 , (21) 



where 



A(zi,ti, z 2 ,t 2 ) = (0\ V {+) (z 1 ,t 1 )^ + \z 2 ,t 2 )\^) = / dhdk 2 (22) 



^ 2y/27lU k2 2y/2lXtU kl ^ u 2 ^ + ^ 1 2 ^ 

Note that the function A(z\, ti, z 2 , t 2 ) satisfies the Klein-Gordon equation 
with respect to Zi,t\ and 22,^2- 

Let us suppose that one of the photons is observed in a frame of reference 
moving with velocity v 1 and another photon is observed in a frame of reference 
moving with velocity v 2 . By using the stationary phase method we obtain 
for large t\ and t 2 : 



A(zi,ti,z 2 ,t 2 ) = J j£- J 2 * J '{k w ,k 20 ) — J: e ™' 2 

hco (ho) V t 2 u (ho) 2y/2iruj kl0 2^2nu k2 



V20 



(23) 

^ e -it 1 (ui kl0 -ik 1Q v 1 ) e -it 2 (u> kQ2 -tk20V2) _|_ ^ ^ fe\ j 
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Here k w = mv\/ \J\ — v\ , k 20 = mv 2 / y/l — v 2 , Z\ = v\t\ and z 2 = v 2 t 2 . 
Therefore 

Ply f y f \ - |/(fcl0,fc 20 ) | 2 , .* 

It is interesting to see the difference between the entangled wave function 
/(&10, ^20) an d the separable one by looking to it with an explicitly indicated 
dependence on the spacetime coordinates: 



z\ z 2 - 2 



f(k 10 ,k 20 ) = /(m^/Wl - ^mfjJl - |) (25) 

Let the wave function of two photons /(fci, k 2 ) in a waveguide be a smooth 
fast decreasing function. Then the probability of observing two photons 
in the waveguide should satisfy the following bounds. For any rii,n 2 = 
0, 1, 2, 3, ... there exist constants C nin2 such that for 

\zi\ > \h\, \z 2 \ > \t 2 \, (26) 

one has 

P(z u h, z 2 , t 2 ) < , l ° ni / n2 ;— p— . (27) 

1 ' ' ' ' " (1+1^1)^(1+ \ Zl \)n 2 ^ I 

Furthermore there exists a constant C such that 

P(z 1 ,t ll z 2 ,t 2 ) < (28) 
|ii| |i2| 

for all zi, ti, z 2 , t 2 . An asymptotic estimate for C, valid for large t\ or t 2 , is 
provided by (1241) . 

An explicit expression for the wave function of the biphotons in a special 
case is given in [2T] : 



f(h,k 2 ) = ^f P {h + k 2 W6hk 2 {h + k 2 ), (29) 
fc 

where fp{k\ + k 2 ) is a Gaussian function describing the pumping photons. 
By using this form of the wave function we obtain the asymptotic formula 
for the probability in this special case. 
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To conclude, the main result of this paper is the bounds to the probability 
density (pQ) and (I2"6"j) . (12"T|) which, in principle, should be possible to test in 
experiments. 
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